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Many results from astrophysical observations point to a 27% contribution of non-baryonic dark 
matter to the mass-energy budget of the universe. Although still elusive, strongly motivated can¬ 
didates in form of weakly interacting massive particles could explain the nature of dark matter, 
and their annihilation or decay would give rise to detectable signatures in gamma-rays. In 2012, 
the H.E.S.S. collaboration started taking data with the largest imaging atmospheric Cherenkov 
telescope in the world which significantly lowered the energy threshold of the already operational 
four-telescope system. In particular, due to its location and improved performance at low ener¬ 
gies, the H.E.S.S. experiment is now in a position to extend the search for dark matter line signals 
down to the 100 GeV mass range. The sensitivity to line searches with a new full likelihood 
approach will be discussed and preliminary results from observations with the second phase of 
H.E.S.S. will be presented. 
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1. Introduction 

Weakly interacting massive particles (WIMPs) are currently among the best candidates to ex¬ 
plain the elusive nature of dark matter (DM) and are actively searched for by the community (see 
[1] for a review). In particular, indirect detection using gamma-rays is considered as one of the 
most promising avenues as it allows the precise determination of the site of their production, and 
high-energy gamma-ray experiments have now the capability to reconstruct accurately the observed 
spectra. The gamma-ray spectrum arising from the annihilation of WIMPs is characterized by two 
main components: a continuum of gamma-rays up to the WIMP mass and possible gamma-ray 
lines. The detection of such a line would be a smoking gun signal for the existence of WIMPs 
since no standard astrophysical process can mimic a line signature. Its search is therefore facili¬ 
tated and much more robust than the search for the continuum, in particular for regions of the sky 
with strong astrophysical gamma-ray emissions [2]. 

The H.E.S.S. collaboration reported on the search for line signatures based on observations of 
the Galactic centre performed with H.E.S.S. I, covering a WIMP mass range between 500 GeV and 
20 TeV [3]. The commissioning of the CT5 telescope lowered the energy threshold down to ~80 
GeV at the analysis level for zenith observations. The new H.E.S.S. II data now allow probing the 
unexplored mass range below 500 GeV and to achieve a significant overlap with the H.E.S.S. I and 
Eermi-EAT results. 

2. Methodology 
2.1 Analysis strategy 

The general analysis strategy relies on the idea of fitting a line-like signal on top of a measured 
background energy spectrum in the ON-source region, using a full event-by-event likelihood pro¬ 
cedure optimized for DM searches in the Galactic Center region. Here no background subtraction 
was performed therefore preserving the entire potential dark matter signal. Since measured distri¬ 
butions are considered in the fit there is no need for acceptance corrections thus strongly limiting 
the associated systematic uncertainties. However, as there is no "universal" background that can be 
estimated from H.E.S.S. observations due to the difference in night sky background distributions 
over a given EoV, a careful choice of background control regions was performed. 

Eor the analysis presented, the region of interest (ROI) is centred at the best fit position of the 
130 GeV line feature reported in Eermi-EAT data ([4],[5]) that is found to be displaced by —1.5° 
longitude with respect to the Galactic Center, called here "Eermi hotspot". The data control re¬ 
gions were defined closely and symmetrically surrounding the Eermi hotspot position to estimate 
the background fraction in the analysis region. Spectral contributions from cosmic-ray background 
and astrophysical diffuse emission are supposed to be included. The instrument response functions 
(IREs) obtained from gamma-ray Monte-Carlo (MC) were used to derive the expected measured 
energy spectra for the line signals considering the stereo reconstruction algorithm which infers the 
direction and energy of gamma-like events from both the signal and the background regions. 
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Sensitivity estimates are provided by MC simulations in section 3 for a line scan between 100 
GeV and 2 TeV, and preliminary results with 95% confidence level (CL) limits are given in section 
4 on a sub-sample of the total H.E.S.S.II dataset focusing on the particular case of the 130 GeV 
line feature. 

2.2 Full event-by-event likelihood 

The full likelihood function is composed of a Poisson normalization term (based on the total 
number of events in the signal and background regions) and a spectral term related to the expected 
spectral contribution of the signal and the background component in the analysis region (ON-source 
region). A description of the method applied in this study is given in [6]. The number of signal 
(Nsignai) and background (Nbckg) events in the analysis region are free parameters of the model, 
while additional information on the signal and background spectral shape is included in the fit. 
Usually, the extracted parameters are the line signal fraction rj = jy which represents the 

relative contribution of the signal in the analysis region, and the line energy position E/i„g. Here 
results on rj are presented. The full likelihood formula is: 


jNsignal+Nbckg)^'^^ „-{N,i„,„l+Nhrk«] V -aNbckg 


(A^ signal j A f,g/cg | A Qbl j A OFF ; A,') 



where Nqn and Noff are the measured number of events in the signal and background regions, a 
the exposure ratio between background and signal regions, and F, represent the individual energies 
of the events in the analysis region. PDFgignai and PDFbckg are the probability density functions for 
the signal and background components that refer to measured energy spectra i.e. photon energies 
smeared by the response functions of the instrument. PDFgignai is a function of the line energy 

2.3 ROI optimization 

An alternative binned likelihood method was developed to perform the ROI optimization. The 
binned ON-OFF excess is fitted to an IRF-convolved power-law in energy side-bands above and 
below the line energy region. The excess due to diffuse emission in the line energy region is 
estimated from the fit. Given the expected and observed number of events in the line energy region, 
the averaged upper limit is computed using Rolke et. al [7]. 

Assuming the peak of the DM density profile at the center of the ROI, the size of the signal region is 
chosen to maximise the sensitivity determined with the binned likelihood method. The optimization 
was done with a 130 GeV line signal following the annihilation cross-section given in [4]. Diffuse 
gamma-ray emission is estimated using the H.E.S.S. Galactic Plane Survey [8], extrapolated to 
lower energies using a power law of index 2.3, while background due to protons and electrons is 
estimated using OFF-data from control samples. An optimal signal region size of 0.4° is found, 
corresponding to a solid angle Afl = 1.531 * lO^^^r. 
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3. Sensitivity Estimates 

The full likelihood method was performed with MC simulations to provide sensitivity esti¬ 
mates as a function of and to study the impact of systematic uncertainties. The numbers of 
measured background events in the ROI of 0.4° and PDF^ckg parametrization were derived from 
the control OFF-source regions around the Fermi hotspot. The form of the PDF^ckg is shown in 
Figure 1 . The PDFsig„ai was provided by a full MC simulation of a line signal at given energy Fum 
smeared by the H.E.S.S. energy resolution. The full likelihood fits covered two pre-defined energy 
ranges from 80 GeV to 1 TeV and from 200 GeV to 3 TeV which allowed probing the sensitivity 
to line signals from 100 GeV to 500 GeV and from 500 GeV to 2 TeV, respectively, ensuring a 
large energy lever-aim in the fit in each case. Sets of 500 simulations were considered for each 
line energy of the scan assuming lOOh of exposure, and led to 95% CL upper limits first on 77 , then 
on the number of excess events The resulting limit is defined as the median of the 500 

simulations, as shown in the example in Figure 2. 



Figure 1: Background PDF for the full likelihood 
method. The PDF is generated from background 
control regions around the Fermi hotspot (see re¬ 
sults in section 4). The line signal PDF at 500 GeV 
is also represented in red for 77 = 0.05. 



Figure 2: Distribution of the 95% CL limit on the 
number of excess gamma-line events in the 

analysis region using the full likelihood method, 
here with the example of a 130 GeV line. The limit 
is taken as the median in the 500 MC datasets. 


The corresponding limits on the flux (<I>) and on the DM annihilation cross-section (< ov >) are 
given by: 

Emax 

^t95%cl f 

<J)95%CL _ y ^ _ Efnin 

ToBS ^max 

f A,ff{Fy)mFy)dEy 

F ^ 

‘-'mm 

where Tqbs is the observation time, Agff and ^ are respectively the effective area for gammas (in 
m^) and the differential energy spectrum of the expected DM signal expressed as functions of the 
true energy, the DM particle mass, [Emin > Emax] the bounds of the energy range and ^astro the 
astrophysical factor. A DM Einasto profile [9] with halo parameters given in [3] has been con¬ 
sidered at the center of the ROI which leads to the value of ^astro = 2.46 x 10^* GeV^cm^^. Eor 
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DM gamma-ray lines the differential energy spectrum is ^{Ey) = 25{Ey — m^) where the factor 
2 results from the annihilation of DM particles into two gammas. The impact of radial acceptance 
within the signal region (not considered in the calculation) has been assessed and is found to only 
affect the limits obtained at the few percent level. The signal region being sufficiently large there 
is no effect due to the point spread function (PSF). Also due to the large extension of the galactic 
DM halo a fraction of the expected DM signal leaks into the background regions, at the level of 
25.1% of the DM signal in the ROT The presented < av > limits account for that effect. Finally, 
nuisance parameters have been introduced in the full likelihood function to estimate the impact of 
systematic uncertainties in the limits calculation. The considered sources of errors are the IRFs, the 
background PDF shape and the diffuse emission component included in the background regions. 

Results on the flux and < av > limits with MC data are presented with the blue data points 
in Figures 3 and 4, respectively, and show the potential of the method for line signal detection. 
The flux is expressed per solid angle unit (with a division by AD) to be compared with the prior 
line analysis at the Galactic Center[3] performed with a larger ROT As an example at 500 GeV, 
^95%cl^^q _ 2 3 X jq-5 ^ y95%CL_ 3 g ^ cm^s^'. The data limits 

(violet) seen on Figure 4 will be discussed later in section 4. The background PDF shape has been 
identified to be the dominant source of systematic uncertainties in the limits calculation, represented 
with a dashed line. The limits obtained with H.E.S.S. II are completing the former H.E.S.S. I [3] 
and Eermi-EAT [10] scans and cover the gap in mass between 300 and 500 GeV. 



Figure 3: Full likelihood results : Flux limits at 95% CL for a line scan between 100 GeV and 2 TeV. The 
flux is expressed in 7 and per steradian. The filled blue data points represent the computed MC 

limits with lOOh of observation time and taking the median value of 500 simulations. Effects on the limits 
due to systematic uncertainties on the background PDF are represented by a dashed blue line. Former limits 
from H.E.S.S. I [3] are represented by red data points. 


5 





Search for Gamma-ray Line Signatures with H.E.S.S. 



Figure 4: Full likelihood results : < dv > limits at 95% CL for the line scan between 100 GeV and 2 
TeV, expressed in cm^s^*. The filled blue data points represent the computed MC limits assuming lOOh of 
observation time and taking the median value of 500 simulations. Effects on the limits due to systematic un¬ 
certainties on the background PDF are represented by a dashed blue line. The data points in violet represent 
limits obtained from 2014 data sample corresponding to 2.8 hours of observation time. For comparison, the 
expected MC limits for the same observation time (rescaled from the 100 hour values) are also shown with 
a violet solid line. Former limits from H.E.S.S. 1 [3] and Fermi-LAT [10] are represented by red and black 
data points, respectively. Finally the < dv > value corresponding the the 130 GeV line feature reported by 
C. Weniger [4] is shown in green. 


4. Particular Case of the 130 GeV Line Feature 

The centre of the 130 GeV excess as observed in the Fermi-LAT data was found to be displaced 
with respect to the position of the Galactic Centre by —1.5° galactic longitude, although with large 
uncertainty [5]. Since the exact position of this excess appears uncertain, H.E.S.S. observations 
were performed in a scan pattern along the Galactic plane, with pointing positions ranging from 
—2.3° <I < 0.5° with a step size of 0.7° and b = ±0.8° in galactic coordinates. 

A total of ~20h of data covering the Fermi hotspot position have been taken in spring and summer 
2014 of which only 2.8 hours have been analysed till now. At least 4 telescopes were requested 
for ensuring good reconstruction of the gamma ray events, including the large CT5 telescope and 
considering observations at low zenith angles (< 20°) to guarantee the lowest possible energy 
threshold. Data quality checks have been performed by checking the global run and the individual 
telescope status. Cuts have been applied on the telescope trigger rates, trigger rate stability and 
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Significance Map Significance distribution 



Preliminary 


Figure 5: Significance map (left) and significance distribution (right) in the FOV. The ROI is represented 
with a white circle centred on the Fermi hotspot (—1.5°, 0°). The known source HESS J1745-290 is clearly 
detected. 


broken pixels fraction in the camera in order to remove bad quality runs from the analysis. Af¬ 
ter data quality selections, several criteria on photon—hadron discrimination were applied to the 
reconstructed event samples in the signal and background regions leading to an effective energy 
threshold of ~ 80 GeV. 

Then significance maps were reconstructed for the 2.8h dataset as shown in Figure 5. The dis¬ 
tribution of significances derived from background fluctuations follow a Gaussian distribution as 
expected, with a significant excess in the Field-of-View at the position of HESS J1745-290. No 
significant excess (< la) is found in the 0.4° ROI at the best-fit position of the Fermi hotspot 
(/,Z>) = (-1.5°,0°). 

Consequently, 95% CL upper limits on the flux and < av > have been derived from the full like¬ 
lihood fit on the data sample following exactly the same procedure as for the MC simulations 
discussed in section 3 and considering the same DM Einasto profile [3]. Limits for the complete 
scan between 100 GeV and 2 TeV are shown in Eigure 4 (violet data points) and are compared with 
the expected limits from the MC simulations, rescaled to 2.8h of observation time. Very good com¬ 
patibility between MC predictions and the data analysis is observed. The related data/MC values 
at 130 GeV are presented in Table E The limit on < av > at 130 GeV exceeds the best fit value 
of [4] value by a factor of 1.3 with ~3 hours of observation time. Only a substantial increase in 
observation time would lead to more stringent limits on the cross-section found by C. Weniger [4]. 
A cross-check study was performed on the source data sample with different calibration and recon¬ 
struction methods. The preliminary results confirm the conclusion of an insignificant excess at 130 
GeV as obtained in the analysis described above. 
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10 ^^ ym^^s^^sr^^ 

< av >95 %cl 

10 ^^^ cm^s^^ 

Data 

9.96 

1.64 

MC 

9.95 

1.64 


Table 1: 95% CL limits on the flux (per solid angle unit) and < (7v > from the likelihood fit for the detection 
of the 130 GeV line. The MC values are coming from the MC line scan presented in section 3, rescaled to 
2.8h of observation time. 

5. Conclusions 

The present analysis of the H.E.S.S. data including CT5 in the setup suggests an excellent 
potential for line reconstruction between ~100 GeV and and a few TeV, considering lOOh of ob¬ 
servation time. The study was performed with background estimated in the Galactic Centre region 
and yielded limits on the gamma-ray flux and on < av > which efficiently fill the gap between 
results from Fermi-LAT and the first phase of H.E.S.S. On the other side, the analysis of a limited 
set of data corresponding to ~3 hours of observation time shows no hint of an excess at 130 GeV 
at 95% CE in the region of the Fermi hotspot {l,b) = (—1.5°,0°). 
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